Rubidium superoxide: a p-electron Mott insulator 
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Rubidium superoxide, Rb02, is a rare example of a solid with partially-filled electronic p states, 
which allows to study the interplay of spin and orbital order and other effects of strong electronic 
correlations in a material that is quite different from the conventional d or / electron systems. Here 
we show, using a combination of density functional theory (DFT) and dynamical mean-field theory, 
that at room temperature RbCh is indeed a paramagnetic Mott insulator. We construct the metal- 
insulator phase diagram as a function of temperature and Hubbard interaction parameters U and J. 
Due to the strong particle-hole asymmetry of the Rb0 2 band-structure, we find strong differences 
compared to a simple semi-elliptical density of states, which is often used to study the multiband 
Hubbard model. In agreement with our previous DFT study, we also find indications for complex 
spin and orbital order at low temperatures. 



I. INTRODUCTION 

Rubidium superoxide, Rb02, is an interesting example 
of a material, where spin and orbital order appears not 
as a result of partially filled d or / states, but due to 
partially filled p electron states. RbC>2 is a member of 
the family of alkali superoxides AO2 (A = K, Rb, or Cs), 
which are insulating crystalline materials composed of 
A + and (02)" ionsi^. At room temperature, RbC>2 has 
a tetragonal crystal structure (see Fig. HJa)), while with 
decreasing temperature this structure undergoes several 
weak distortions, first to orthorhombic, then to mono- 
clinic symmetry 2 ^ 3 -. The electronic structure around the 
Fermi level is dominated by oxygen p states which can 
be well approximated by molecular orbitals (MOs) corre- 
sponding to the O2 units, and are filled with 9 electrons 
(see Fig.[ljb)). Assuming no further symmetry breaking, 
the two highest occupied antibonding ir* orbitals are 3/4 
filled. 

The degeneracy of these orbitals can be lifted through 
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FIG. 1: (a) Tetragonal crystal structure of Rb02 at room 
temperature, (b) Electronic structure represented by a set of 
oxygen p molecular orbitals. 



cither magnetic or orbital long range order, or both. The 
alkali superoxides thus allow to study "correlation ef- 
fects" in a completely different class of materials com- 
pared to the more conventional transition metal oxides 
or / electron systems. Antiferromagnetic order is indeed 
found experimentally at low temperatures (Tn(Rd02) « 
15 K)i*2, and it was suggested by recent density func- 
tional theory (DFT) and model studies that the insulat- 
ing character of alkali-supcroxides at low temperatures 
can be explained by the interplay of correlation effects 
(spin and orbital order) and crystal distortions^—. How- 
ever, the nature of the insulating state of these superox- 
ides at room temperature has so far remained unexplored. 

Due to the high symmetry crystal structure with no 
long-range order of spins or orbitals, it is impossible to 
explain the insulating character of the alkali superoxides 
at room temperature within an effective single particle 
band picture. Here we show, using a combination of DFT 
and dynamical mean field theory (DFT+DMFT), that 
Rb02 at room temperature is in fact a Mott insulator, 
where the strong Coulomb repulsion prevents the electron 
hopping between adjacent sites. 



II. ELECTRONIC STRUCTURE OF RB0 2 

We obtain the electronic structure of Rb02 from a 
non-spinpolarized DFT calculation using the Quantum- 
ESPRESSO package^, employing the generalized gra- 
dient approximation of Perdew, Burke, and Ernzerhofii 
and ultrasoft pscudopotentialsi 2 -. Figs.^a) and (b) show 
the resulting density of states (DOS) and bandstructure. 
It can be seen that the electronic structure of Rb02 in- 
deed closely resembles the simple MO picture sketched 
in Fig. QJb), with a splitting of about 5eV between the 
bonding and antibonding 7r and 7r* bands, and a sin- 
gle band corresponding to bonding a MOs at — 6eV. 
The antibonding a* states at ~ 5 cV arc strongly in- 
termixed with other empty states corresponding to the 
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FIG. 2: (Color online) Density of states (DOS) (a) and band- 
structure (b) of nonmagnetic Rb02. In (a) the total DOS 
is shown as thick (black) line, while the projection on O2 a 
and 7r states are shown as (blue) striped and (red) filled ar- 
eas, respectively. Panel (c) shows the DOS (per spin-orbital) 
of the 7T* bands used in the DMFT calculations, compared 
to a semi-circular DOS with the same bandwidth (a small 
broadening is applied in both cases). 



Rb + cations. 



III. DYNAMICAL MEAN-FIELD THEORY — 
COMPUTATIONAL METHOD 

To calculate the electronic properties at finite temper- 
ature and account for local correlation effects, wc use 
dynamical mean- field theory (DMFT) which allows to 
map the lattice problem to an effective problem of a 
single-site impurity surrounded by a bath. The inter- 
action part of the impurity Hamiltonian is taken to be of 
the Slater-Kanamori form 

~ J ( d l,l d lt db >-t da ^ + 4,t4,4. d «.t d «4 + h - c -)> (!) 

a=/tb 

with d\ a the creation operator for an electron of spin 
a in orbital a and U' = U — 2J. To solve the ef- 
fective impurity problem, we use the strong-coupling 
continuous time quantum Monte Carlo approach (CT- 
HYB)— liil. From the self-consistcntly determined hy- 
bridization function A(t), the impurity Green's function 
Gi mp (r) is computed and measured on a homogeneous 

^//3/40eV _1 



grid of N T = 1000- 



points, where [. . . ] rep- 



formation we obtain the self-energy in Matsubara space, 

S(iw„) = iuj n +n - Gfa p (iw n ) - A(iw„), (2) 

where ui n = (2n + 1)tt//3 for integer n, /i is the chemical 
potential, and (3 = 1/T the inverse temperature. Using 
this self-energy and the single-particle Hamiltonian H(k) 
we obtain the local lattice Green's function by averaging 
over the Brillouin zone: 



(iw n ) = — [iuj n +fj,- H(k) - E(iw n )] 1 . (3) 



The DMFT self-consistency condition demands that this 
local lattice Green's function is the same as the impu- 
rity Green's function. This condition, in combination 
with Eq. yields the hybridization function for the next 
DMFT iteration, 



A(iw„) 



+ fi- G lo J(iw„) - £(iw„). 



(4) 



resents the nearest integer number. After Fourier trans- 



We only include the partially-filled antibonding ^*/ y 
bands in our DMFT calculations for Rb02, and ex- 
press the corresponding Hamiltonian H(k) in a basis 
of maximally localized Wannier functions^. The cor- 
responding DOS is shown in Fig. [U[c). One can recog- 
nize a pronounced asymmetry with respect to half-filling. 
In DMFT studies, a model semicircle (SC) density of 
states (DOS) is often employed to represent the electronic 
bands, since it leads to a simple expression connecting A 
and G; mp . Furthermore, due to the resulting particle-hole 
symmetry, only occupations between zero and half-filling 
need to be studied. Here, we investigate the differences 
between results obtained using the model SC DOS and 
the realistic DFT band structure of Rb02 in the tetrag- 
onal crystal structure. The bandwidth of the SC DOS is 
set equal to the bandwidth of the Rb0 2 DOS (0.93 eV, 
see Fig. Etc)). 



IV. RESULTS 
A. Room temperature properties 

Generally, the value of the spectral function at zero en- 
ergy indicates whether a material is insulating or metal- 
lic. However, obtaining the spectral function from the 
imaginary time Green's function requires an analytic con- 
tinuation to the real axis, which can introduce additional 
uncertainties^. We therefore consider several possible 
indicators for the metal-insulator transition (MIT) which 
are directly accessible from G; mp (r) . All of these quan- 
tities (described in more detail below) are compared in 
Fig. [31 calculated at room temperature (f3 = 40 eV -1 ) 
for interaction parameters close to the MIT (U = 1.2 eV 
and J = 0). 

One possibility is to monitor the occupation 
n = —^2 a Gf {P) (a is the spin-orbital index) as a 
function of the chemical potential \x (see Fig. EJa)) and 
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FIG. 3: Various quantities evaluated from the impurity 
Green's function indicating the metal-insulator transition at 
integer filling, calculated for /3 = 40 eV -1 (T w 290 K), 
U = 1.2 eV, and J = 0.0 eV. The average over all spin-orbitals 
is shown in panels (b), (c), and (d). 



to identify the insulating phase by a plateau in 
This, however, requires a large number of calculations 
for slightly different values of (i. Another indicator is 
given by the mass enhancement in the low-temperature 
metallic phase, which grows rapidly as the Mott insu- 
lating state is approached, and which we estimate from 
the self-energy at the lowest Matsubara frequency (see 
Ref.01): 
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(5) 



We also consider the following estimate of the spectral 
function (see e.g. Ref. [l8l): 



A(0) 



P Q. 



(6) 



While for both of these quantities (shown in Fig. EJb) 
andGOc), respectively) only one calculation at the correct 
fi value is required, the identification of the MIT phase 
boundary requires the definition of a suitable thresh- 
old value. In addition, G(/3/2) suffers from significant 
statistical noise in the insulating state, as this r-region 
is difficult to sample with standard CT-HYB. A quan- 
tity which is quite insensitive to noise is the slope of 
G; mp (iw — > 0), which is positive/negative for the metal- 
lic/insulating stated. In practice, we estimate the slope 
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FIG. 4: Slope K for Rb02 input and selected values of U and 
J at P — 40 eV -1 (average over all spin-orbitals is shown). 



K from Gi mp (iw) at the two lowest Matsubara frequen- 
cies: 
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(7) 

In Fig. |31 a clear difference can be seen between the 
SC and the RbC>2 input. For the SC, the insulating state 
already appears at 1 and 3 el filling, while 2 el filling is 
clearly metallic. In contrast, the RbC>2 electronic struc- 
ture yields a clear insulating state at 3 el filling and a 
(just barely) insulating state at 2 el filling, while 1 el 
filling is still metallic. The different range of /i which 
leads to the insulating state for SC and Rb02, respec- 
tively, indicates a sizable shift of the corresponding MIT 
boundary. The particle-hole asymmetry in the real elec- 
tronic structure of Rb02 thus leads to large quantitative 
changes compared to the simple SC DOS. 

Fig. 2] shows the slope K for different values of U and 
J at T w 290 K using the Rb02 band structure. For 
fixed J = eV, there is an obvious tendency towards the 
insulating state with increasing U for all integer fillings, 
as expected. However, at fixed U = 2 eV, increasing J fa- 
vors the Mott insulator at half-filling (2 el) but favors the 
metallic solution for 1 and 3 el fillings. This is consistent 
with previous discussions of multi-orbital model a 19 i 20 : In 
the large- U limit, the width of the "Mott plateau" in ^ is 
given by A^ ott = E n+1 + £„-i - 2£„, with E n denoting 
the lowest eigenvalue of the n-particle eigenstates of H mt 
(Eq. (HJ))^. In our two-orbital case this estimate yields 
U — 3J for n = 1, 3 and U + J for n = 2, in agreement 
with the observed dependence of the plateau- width on J. 
In reality, the Mott plateau will be reduced by approx- 
imately the bandwidth W, so that we obtain the rough 
estimate 



A Mott 
^3 



U-3J-W . (8) 
Based on the identification of the MIT boundary us- 
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FIG. 5: Metal-insulator phase diagrams, (a) Critical U as 
function of J at room temperature (T ~ 290 K). The green 
star indicates the realistic values for U and J calculated irA 
(b) Critical temperature as function of U — 3 J. 



ing the slope K, we computed the phase diagram for 
3 el filling at room temperature (T « 290 K) as a 
function of the interaction parameters U and J (sec 
Fig. [5{a)). Since AMott must be larger than zero for a 
Mott insulating solution to exist, Eq. © also provides 
a crude estimate for the critical interaction strength: 
U c — 3J « const. It can be seen that the MIT bound- 
ary in Fig. EJa) agrees nicely with this simple estimate. 
Furthermore, the critical U — 3J at room temperature 
differs by « 30% of the bandwidth between Rb02 and 
the simple SC DOS, and increases slightly as a function 
of T for T > 300 - 400 K, while the opposite trend is 
observed at lower T (Fig. [5{b)). 

For realistic values of the interaction parameters U = 
3.55 eV and J = 0.62 eV, which were obtained for it* 
orbitals in the very similar material KO2 using the con- 
strained LDA and random-phase approximation^, the in- 
sulating state is obtained for both SC and Rb02- Our 
results therefore predict a Mott insulating state (without 
long-range order) for Rb02 at room temperature, consis- 
tent with experimental observations. 

While we did not find evidence for a coexistence re- 
gion, indicative of a first order MIT, for T > 145 K, 
we have verified that we obtain an insulating state with 
a clear gap in the spectral function, even at room tem- 
perature. To demonstrate this, we have used the max- 
imum entropy method^ to construct spectral functions 
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FIG. 6: Spectral functions (a) and the various indicators 
for the MIT (b) for Rb0 2 at /3 = 40 eV" 1 (T « 290 K), 
J = 0.6 eV, and different values of U around the MIT. 
Z = (m*/m)~J and K in eV -2 . The noninteracting DOS 
is shown as grey shaded area in (a). 



at p = 40eV _1 , J = 0.6 eV, and different values of U 
around the MIT. The result is shown in Fig.^a). A gap 
is present for U > 2.7 eV, in perfect agreement with the 
various indicators of the MIT discussed previously (and 
which are shown in Fig.H^b)). In agreement with Ref. l2ll 
we find a "bad metal" region with a strongly renormal- 
ized Z — (m* /m)^ 1 < 0.4 in the vicinity of the Mott 
transition. From Fig. EJa) it can be seen that this cor- 
responds to spectral functions with substantial narrow- 
ing of the central quasiparticlc feature and an emerging 
three-peak structure visible for U = 2.4 eV. In addition, 
there is a significant spectral weight transfer to energies 
around — 2eV compared to the noninteracting DOS. 



B. Low temperature behavior 

Finally, we focus on the low-temperature behavior. 
While for temperatures T > 200 K (for which wc didn't 
find indications of ordered states), the hybridization 
function is averaged over all spin-orbitals in each iter- 
ation, Fig. [7] shows the evolution of the occupation of 
each individual spin-orbital at T 29 K when no such 
averaging is performed. While for some values of J the 
occupation eventually converges to spin and/or orbitally 
polarized states, the occupations exhibits characteristic 
oscillations for other values of J. As discussed ir>22, such 
oscillations indicate that the system wants to adopt an 
ordered state with a sublatticc structure that is incom- 
patible with the applied self-consistency condition (in our 
case all sites are forced to be equivalent). 

Even though we do not attempt to fully resolve the re- 
sulting spin and orbital patterns, we can make a number 
of interesting observations. First of all, there are dras- 
tic differences between Rb02 (Fig. [7{a)) and the simple 
SC DOS (Fig. [Tib)). The latter oscillates between three 
states with different spin and orbital polarization (SP and 
OP) and is insulating for J < 0.8 eV, while for higher J 
it is metallic with no SP and OP (in these calculations 
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and almost no OP. For J = 0.62 and 0.8 eV, a stable 
solution with large OP and zero SP appears, while a fur- 
ther increase of J induces a stable SP and reduced OP. 
The system is insulating for all J < 1.0 eV, while for 
J = 1.2 eV it is a ferromagnetic half-metal with full 
SP and no OP. In Fig. [gja) we show the SP and OP of 
Rb02 as a function of temperature for the realistic values 
U = 3.55 eV and J = 0.62 eV. The system is insulating 
and while essentially no SP develops down to T ps 30 K, 
OP appears below T m 60 K and reaches almost its max- 
imum at T ps 30 K. While it is not possible from our 
calculations to make a prediction about the character of 
the expected spin- and orbitally-ordered ground state, 
the above temperatures are consistent with our previ- 
ous estimate of the ordering temperature based on total 
energy differences of different orbitally ordered configu- 
rations obtained from DFT+C/ calculations at T — K— . 
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FIG. 7: (Color online) Evolution of individual occupations for 
different J and fixed U - J = 2.93 eV at T fa 29 K for Rb0 2 
input (a) and for the SC DOS (b). (Black) Squares: orbital 
1, spin up; (red) diamonds: orbital 2, spin up; (green) plus 
symbols: orbital 1, spin down; (blue) crosses: orbital 2, spin 
down. 
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FIG. 8: (a) Temperature dependence of spin and orbital po- 
larization for fixed U and J. (b) spin and orbital polarization 
as a function of J for U - J = 2.93 eV and T ps 29 K. All 
data corresponds to Rb02 input. 

both U and J have been varied while keeping U c s = U—J 
constant). For Rb02 we can distinguish three different 
regimes (see Fig. [5£b)). For J < 0.5 eV, the occupation 
oscillates between three different states with different SP 



V. SUMMARY AND CONCLUSIONS 

In summary our calculations clearly show that for re- 
alistic values of the interaction parameters U and J, 
Rb02 at room temperature is a paramagnetic Mott in- 
sulator without exhibiting any symmetry-breaking long- 
range order. We find pronounced quantitative differences 
between the widely used SC DOS and the realistic elec- 
tronic structure of Rb02 , which leads to a strong asym- 
metry between the 1/4-filled and 3/4- filled cases. We 
also find indications of complex spin and orbital order 
below T ps 30 K, the character of which seems to depend 
strongly on J. Furthermore, at low temperature Rb02 
exhibits clear qualitative differences compared to the sim- 
plified SC DOS. It will be interesting to clarify in future 
work whether single site DMFT is capable to resolve the 
complicated spin and orbital patterns predicted within 
model calculations based on a perturbative treatment of 
electron-electron interaction and a simplified electronic 
structure of RbO^i^- 
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